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ABSTRACT 

We calculate the location of the magnetically-inactive dead zone in the 
minimum-mass protosolar disk, under ionization scenarios including stellar X- 
rays, long- or short-lived radionuclide decay, and energetic protons arriving from 
the general interstellar medium, from a nearby supernova explosion, from the disk 
corona, or from the corona of the young star. The disk contains a dead zone in all 
scenarios except those with small dust grains removed and a fraction of the short- 
lived radionuclides remaining in the gas. All the cases without exception have 
an "undead zone" where intermediate resistivities prevent magneto-rotational 
turbulence while allowing shear-generated large-scale magnetic fields. The mass 
column in the undead zone is typically greater than the column in the turbulent 
surface layers. The results support the idea that the dead and undead zones are 
robust consequences of cold, dusty gas with mass columns exceeding 1000 g cm -2 . 

Subject headings: circumstellar matter — solar system: formation — stars: for- 
mation — instabilities — MHD 



1. INTRODUCTION 



Angular momentum transport is a key to the evolution of protostellar disks and the 
origins of the planets, as it governs the flow of raw materials toward the star (IBodenheimer 
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19951 ). Transport processes likely to be important include the turbu lence resulting from the 



magneto-rotational instability (MRI; iBalbus fc Hawlev 



large -scale magnetic fields driving an outflow ([B 



19931 ) or shearing within the disk ([Turner fc Sano 



1991 



1998) and the stresses due to 



20081 ) 



andford fc Payne! Il982l ; I War die fc Konigl 



All of the magnetic angular momentum transport processes require the gas is sufficiently 
ionized to couple to the fields. Collisional ionization leads to good coupling only at the tem- 
peratures above 1000 K found very close to the star (IPneuman fc Mitchelllll965l ; lUmebayashi 
19831 ; iNakano fc Umebayashilll988l ). In the overwhelming majority of the disk volume, the 
main ionization processes are non-thermal . Among the import ant non-thermal processes 
are the X-rays emitted by the young star (|Glassgold et al.lll997D. the c osmic rays arriving 



from interstellar space (lUmebavashi 



1983 



Umebavashi Sz Nakanolll988l ). and the decay of 



radionuclides within the gas ( Umebavashi &: Nakanol 1981 ) . The basic difficulties in reachin. 



high enough levels of ionization through these non-thermal processes are the large column of 
absorbing material and the high rate of recombination on the surfaces of dust particles. The 
X-rays and cosmic rays generally penetrate and ionize only the surface lay ers of the disk , 



leaving a region near the midplane whe re magnetic activity is suppressed (iGammid 11996 
Sano et aliEoOol : llkner & Nelsonll2006ah . 



The absorption of X-rays and cosmic rays in the atmosphere means that the disk can 
be divided into three zones. In the upper layers, and in outer annuli where the mass column 
is low enough for the ionizing radiation to reach the midplane, the magnetic fields couple 
thoroughly to the gas, and magnetic forces drive turbulence through the MRI. In a "dead 
zone" near the midplane, the fields decouple from the gas, and magnetic forces are largely 
irrelevant. At resistivities between these extremes, magnetic fields decouple over scales com- 
parable to the disk thickness, shutting off the MRI, but remain coupled over the disk radius. 
The radial gradient in orbital frequency can then shear out any weak radial magnetic field to 
gener ate toroidal fields, enabling angular momentum transport to continue (ITurner fc Sano 
20081 ). The turbulent surface layers are thus separated from the midplane dead zone by an 
"undead zone" that becomes magnetically active when supplied with radial fields. 

In this paper we explore whether three less well-studied additional sources of ionization 
can reduce the sizes of the dead and undead zones or make protostellar disks magnetically 
active throughout. T he three ionization s ources are high-energy particles from (1) a nearby 
supernova explosion (IFatuzzo et ah] 120061 ). (2) the corona of the protostellar disk itself, and 
(3) the young star. We estimate the ionization rates to order of magnitude (£J2]) and calculate 
the resulting resistivities (§3]) in the minimum-mass model of the protosolar disk (§!]), finding 
the undead and dead zones shown in $51 A summary and conclusions are in §6j 
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2. IONIZATION 

The base ionization rate in our calculations is set by the X-rays emitted from the vicinity 
of the young star ( §2.ip . Some of the models include also the decay of radioactive isotopes 
within the disk ( §2.21) . Our purpose is to explore the additional effects of energetic protons 
from three sources: the interstellar medium ( §2.31) . the disk corona ( §2.41) . and the stellar 
corona ( §2.5)) . With the aim of determining whether the dead zone can be removed under 
favorable circumstances, we generally err on the high side in estimating the proton fluxes. 
The rates of ionization from the five processes are compared in §2.61 



2.1. Stellar X-Rays 



The X-ray ionization is included because a lmost all lightly-obscured young solar-mas s 
stars in the Orion nebula show X-ray emission (IGarmire et al.ll2000l ; iPreibisch et al.ll2005l ). 



The measured spread in luminosity is several decades, with a median around 2 x 10 erg s 
Evidence that the X-rays strike the disks comes from the iron Ka flu orescent emission de 



tected in some pre-main-sequence stars having near-infrared ex cesses (jlmanishi et al.l 12001 



Favata et al.ll2005al ; iTsujimoto et al.ll2005l ; iGiardino et al.l 120071 ). Gas at distance r directly 
facing the star is ionized at a rate (x = vFx/Ei « 5 x 10 -11 s _1 (r/AU) -2 , using a flux 
Fx corresponding to the median X-ray luminosity, an energy requirement per ionization 
Ei = 37 eV, and an ab sorption cross-section a = 4.16 x 10~ 24 cm 2 at photon energy 5 keV 
and Solar composition (IGlassgold et al.lll997l ). Most of the disk however is ionized at much 
lower rates, because the photons enter at grazing angles and only the most energetic and 
least numerous reach the deep in terior. Detailed Monte Carlo radiative transfer results for a 
5 keV thermal source spectrum (jlgea fc Glassgoldlll999l ). fitted assuming an inverse square 
falloff with the radius, yield an ionization rate 



(x = 2.6 x 10- 15 s- 1 (r/AU)- 2 [exp(-E a /E x ) + exp(-E & /E x )] 



(1) 



at columns greater than 1 g cm 2 , where Ex = 8.0 g cm 2 is the X-ray absorption depth 
and E a and Ej are the mass columns lying vertically above and below the point of interest. 



2.2. Radionuclides 



The decay of long-lived radionuclides, primarily 40 K (half-life 1. 25 Gyr), gives an ioniza- 
tion rate (lr — 6.9 x 10~ 23 s _1 at protosolar elemental abundances (jUmebayashi fc Nakano 



198ll ). Short-lived radionuclides, primarily 26 Al (half-life 0.717 Myr), yield an ionization rate 
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Csr = 3.7 x 10~ 19 s -1 ( 1Stepinskilll992T) . The results of our cal culations change little if we use 
instead the rates advocated by lUmebayashi &: Nakand (120091 ). which differ from these by be- 
tween 9% and 55%. The radionuclide ionization rates will be substantially lower if the dust 
particles containing most of the heavy elements are absent. In some cases without grains 
we include no radionuclide ionization, while in others the rate is reduced 10 000-fold from 
the dusty calculations to account for the possibility that a small amount of the radioactive 
element occurs in the gas phase. 



2.3. Interstellar Cosmic Rays 



The first of three energetic proton sources we consider is the interstellar cosmic rays. 
Particles with energies above .1 GeV provide most of the ionizing effect at high mass columns 
( lUmebayashi fc Nakanolll98ll ). Shielding by the contemporary Solar wind reduces the flux 
of interstellar cos mic rays reaching the top of the Earth's atmosphere in the energy range 
up to a few GeV ( Spitzer fc Tomaskol 1968h. Estimates of the cosmic ray energy spectrum 
ou tside the heliosphere, re viewed by llndriolo et al.l (120091 ) , scatter around the levels inferred 
by IComstock et al.l (Il972l ) that we use here. The intensity 0.94(Eq + E)~ 2S cm~ 2 s _1 sr _1 
GeV -1 falls off with energy E above the proton rest energy E = 0.938 GeV. The particles 
with E > 0.1 GeV deliver an energy flux of approximately 0.01 erg cm -2 s _1 , producing an 
ionization rate 



-18 -1 



(cr = 5 x 10- 18 s 



exp 



J CR 



J CR 



+ exp 



^CR 



J CR 



(2) 

under isotropic illumination of the disk surface (lUmebayashi fc Nakand 120091 ). The cosmic 
ray absorption depth Y>cr = 96 g cm -2 ( lUmebayashi fc Nakand 1 19811 ). Lacking good infor- 
mation about the degree to which the protosolar wind screened the inner protosolar disk from 
interstellar cosmic rays, and seeking a lower bound on the size of the undead and dead zones, 
we assume that the cosmic rays above 0.1 GeV reach the disk surface with the interstellar 
flux and energy spectrum. 

A supernova explosion occurring within 10 pc of the young star could increase the 
cosmic ray energy density by about three orders of magnitude over the general interstellar 
value, based on the en ergetic-particle popu lations inferred from gamma-ray observations 
of supernova remnants (IFatuzzo et al.l 120061 ) . We therefore consider below the effects of a 
thousand-fold increase in the cosmic ray ionization rate. 
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2.4. Disk Corona 



The second energetic proton source we consider is the atmosphere of the accretion disk 
itself. Forbidden optical line emission and central Balmer absorptions in the spectra of 
T Tauri stars indicate the presenc e of low-density overlying gas with temperatures in the 
range 5 000 to 10 000 K fcwanll997h similar to the Solar chromosphere. Models of the heating 
and cooling balance show that such temperatures can be produced by the stellar ultravio 



let a nd X-ray irradiation (IKamp fc Dullemondl |2004| ; iGlassgold et al.l 12004 ; iNomura et al. 



2001 


; IPvo et al. 


2002 





surrounding a faster central jet. Since magnetic forces likely drive the winds, we assume 
below that the Alfven speeds approach the escape speed in the extended atmospheres of the 
disks. 

Under similar conditions in the Solar chromosphere and corona, fast magnetic recon- 
nection and shock s lead to the acceleration of protons to high energies, as reviewed by 
Miller et al.l (119971 ) . Fast reconnection requires pla sma in the Hall regim e, so that the whistler 
or kinetic Alfven driving waves are not dissipated (jCassak et al.ll2005l ). The electron- ion col- 
lision frequency v e i must be less than the electron cyclotron frequency Q ce . In the atmosphere 
of the disk around a young Solar-mass star, the ratio of the two frequencies is 



a 



0.01 



n 



10 10 cm" 3 



T 



10 4 K 



-3/2 



1AU 



1/2 



(3) 



using expressions from iMitchner fc Krugerl (119731 ) with the Coulomb logarithm set to 10. 
Fast reconnection can thus occur in the tenuous disk atmosphere. Reconnection in the 
dense, cold disk interior presumably is due to the slower Sweet-Parker mechanism. 

Reconnection in the Solar corona produces a mix of thermal and non-thermal particles. 
The temperature of the thermal component is such that the downstream gas pressure is 
approximately equal to the upstream magnetic pressure. Similar conversion of magnetic 
energy to heat in the disk corona yields the virial temperature 10 eV (r/AU)^ 1 . Reconnection 
can thus pla usibly contri bute to the heating and stirring indicated by the central Balmer 
absorptions (IKwanl 1 199 71) and by the Doppl er-broadened emission line of ionized neon at 
12.8 /im in TW Hydrae (IHerczeg et al.l 120071 ). The thermal particles do not affect the dead 
zone as they have little penetrating power. 

The maximum energy of the non-thermal particles is approximately equal to the product 
of the particle charge q, the electric field strength, and the length over which the particles 
are accelerated. The electric field |v x B|/c is proportional to the reconnection speed |v|, 
which we set equal to the Alfven speed Va- The lengths of magnetic loops in the disk corona 
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depend on the rate of reconnection relative to shear (jUzdensky fc Goodmanll2008l ). When 
shear dominates, loops may reach sizes comparable to the radius r, which we choose as the 
acceleration distance. The maximum energy of the particles produced during reconnection 



100 GeV 



n 



10 io cm -3 



1/2 



(4) 



so the high-energy tail extends into the regime of penetrating cosmic rays above 0.1 GeV. 

The luminosity of the particles accelerated by reconnection in the disk corona is bounded 
above by the rate of gravitational energy release due to accretion. The accretion energy per 
unit area and unit time leaving each side of the disk around a young solar analog T Tauri 
star is about 3 x lO 5 (M/lO- 8 M yr- 1 )(r/AU)" 3 times the energy flux in ordinary interstellar 
cosmic rays, where M is the mass flow rate. Much of the accretion flux takes the form 
of thermal radiation from the disk photosphere, while some fraction is in magnetic fields 
that rise into the cor ona before dissipating. Results from shearing-box MHD calculations 
(IMiller fc Stonell2000l ) suggest the fraction is of order 10%. Some part of the coronal magnetic 
energy is converted through reconnection and shocks into the kinetic energy of non-thermal 
protons. We assume the proportion is similar to the contemporary Sun, where magnetic 
energy is dissipated at approximately 10 27 erg s _1 , judging from the mean X-ray luminosity. 
The Sun and solar wind produce a flux at 1 AU of about one proton with energy exceeding 
0.1 GeV, per square centimeter per second, comparable to the interstellar cosmic ray flux 
outside the heliospher e, according to data from the GOES satellite energetic particle sensors 



(IMewaldt et a 



200 7|) and to analy ses of tracks produced in lunar rocks over the past few 



million years (IReedy &: Martil Il99ll ) . The corresponding solar luminosity in particles above 



0.1 GeV is around 10 24 erg s _1 , or 0.1% of the rate of dissipation of magnetic energy. We 
thus infer that the overall fraction of the disk accretion energy that rises into the corona in 
the form of magnetic fields and is converted to energetic protons is 10~ 4 . Perhaps half the 
particles travel downward and strike the disk surface, yielding an ionisation rate Cdp that is 
10(r/AU) -3 times the demodulated interstellar cosmic ray value. 



2.5. Stellar Corona 



The third source of energetic protons we consider is the corona and inner wind of the 
young Sun. The contemporary Sun accelerates ions to cosmic ray energies. The reconnection 
of opposing magnetic fields with Alfven speed va produces ions of mass m with energ ies of 
order mv\, according to o bservations near 1 AU in the So l ar wind (IGosling et al.ll2005f) and 



to particle-in-cell models (IKrauss-Varban fc Welschl 120071 ; iDrake et al. 



gions in the Solar corona have densities of 10 9 to 10 10 cm 3 (iBenz et al. 



1983 : 



2009a ). Flaring re- 



Krucker et al. 
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20081 ) and fields around 50 G, yielding protons with energies in the tens of keV range. 
How these seed particles reach much higher energies remains unclear. Ideas include re- 



peated shock crossings (as reviewed by iMa" 



flLazarian fc Opherll2009l : brake et al.ll2009tt ). 



kov &: Druryl l200ll ) and magnetic reconnection 



The Sun produced more energetic particles when young: the X-ray flares observed in 
young stars, scaled by the solar relationship between X-ray luminosity and proton emission, 
i ndicate a 10 5 -fold enh ancement in energetic protons compared to contemporary solar levels 
( jFeigelson et al.l 120021 ). Such an enhancement appears to be consistent with a spallation 
origin for 26 Al and ot her short-lived radion uclides whose decay products occur in primitive 
chondritic meteorites (jGounelle et al.ll2006l ). In addition, the magnetic fields of young Solar- 
mass stars are strong enough for reconnection to produce large numbers of seed particles with 
energies above an MeV. T Tauri stars have kilogauss s urface-averaged magneti c fields, three 



order s of magnitude stronger than the modern Sun (jJohns-Krull et al.l 119991 ; iJohns-Krull 
20071 ). while the densities of gas on coronal mag n etic loops, inferred from the emitted X- 
rays, are around 10 10 cm -3 (IFavata et al.ll2005bl ; iJardine et al.l 120061 ). Absorption of the 



X-rays by foreground coro nal material indicates columns also consistent with such densities 
(IRobrade fc Schmittl 120071 ). The characteristic energy of the seed protons produced in re- 
connection under these conditions is 10 MeV if the fields are uniformly-distributed, or larger 
if as is likely the fields are patchy across the face of the star. 

Lacking better knowledge of the arrangement of the magnetic fields near the star, we 
cannot say whether the stellar coronal protons are channeled back into the star, out along 
open field lines into a wind, or toward the disk which may reach or penetrate the boundary 
of the stellar magnetosphere. Since our objective is a minimum size for the dead zone, we 
shall assume that the protons propagate in straight lines like the X-rays, reaching all parts 
of the disk surface unimpeded. 

One must also bear in mind that most of the protons are produced in flares lasting 
only hours or days. The flare particles will reduce the size of the dead zone in the disk if 
the ionization created in one flare persists till the next. Calculations of the ionization due 
to X-ray flares suggest that the recombination time is less than th e interval between flare s 
unless the abundance of the dust in the disk is greatly reduced (lllgner &: Nelson! 12006a ). 
A more conservative estimate of the contribution of the stellar energetic particles comes 
from scaling the modern Sun up to the X-ray luminosity corresponding to (x- The effective 
enhancement in the energetic protons is then only about 3.3 decades and the rate of ionization 
is comparable to that resulting from the stellar X-rays. 

Nevertheless, in the interest of obtaining a minimum dead zone size, we treat the 10 5 -fold 
enhancement in the energetic protons as time-steady. The flux falls off roughly as the inverse 
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square of the distance, at distances much greater than the size of the stellar magnetosphere. 
Projection effects reduce the flux per unit disk area by about a factor ten, since the disk 
surface is inclined about 0.1 radian from the line of sight to the star. We take a surface 
ionisation rate (sp of 10 4 (r/AU)~ 2 times the interstellar cosmic ray value. The particles 
from the coronae of the disk and star are assumed to decline with column in the same way 
as the interstellar cosmic rays, following eq. [2l 



2.6. Ionization Processes Compared 

The relative strengths of the different ionization processes are shown in figure [T] where 
the ionization rates at a column of 8 g cm -2 are plotted versus radius. Compared with the 
stellar X-rays, the stellar energetic protons are 40 times stronger and are the strongest process 
plotted. The disk coronal energetic protons are a few times weaker than the X-rays at 0.1 AU 
and fall off faster with distance. The interstellar cosmic rays are distance-independent and 
are stronger than the X-rays outside 16 AU. If enhanced a thousand times due to a nearby 
supernova, interstellar cosmic rays are stronger than the stellar X-rays outside 0.5 AU and 
are the strongest of all the processes outside 3 AU. 



RESISTIVITY 



Recombination in the gas phase and on grain surfaces balances the ionization processes 
in local chemical equ ilibrium. We find the ba lance by solving the reduced chemical reaction 
network discussed by lllgner &: Nelson! (j2006al ). The ionization produces molecular ions which 
are destroyed through dissociative recombination, or by charge exchange with metal atoms 
which then recombine radiatively. Magnesium is chosen as a representative gas-phase metal 
because of its combination of high abundance and low desorption temperature. The amount 
initially in the gas is chosen to be 1% of the Solar magnesium abundance. Recombination 
also occurs on grains which can collide with electrons, ions and other grains. The grains are 
1 fim in radius with an internal density 3 g cm -3 and a mass fraction 1% relative to the 
gas. In addition, molecules and atoms can become adsorbed on grain surfaces, reducing the 
gas phase abundances. The stiff set of ODEs representing the chemical kine tic equations is 
integrated to equilibrium using semi-implicit extrapolation (jPress et al. 1992 ). The resulting 
electrical resistivity rj = 234-v/T / x e cm 2 s _1 varies inversely with the electron fraction x e = 
n P Jn n , where n P is the electron number density and n n the total number density of neutrals 
jBlaes fc Balbu3ll994h . 
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Fig. 1. — Ionization rates at column 8 g cm -2 versus the distance from the star. The dashed 
line shows the contribution of stellar X-rays, and the dotted lines the decay of short-lived (upper) 
and long-lived (lower) radionuclides. Among the three solid lines, the uppermost is for stellar 
energetic particles, the steepest shows the disk coronal particles, and the horizontal solid line is 
for interstellar cosmic rays. Large surface density is assumed and only radiation striking the top 
of the disk is included on this diagram. Circles indicate the results of the Monte Carlo X-ray 
transfer calculations by llgea & Glassgold jl999l ). scaled to a stellar X-ray luminosity 2 x 10 



30 



erg s . The disk coronal particle estimate corresponds to a fraction 10 of the energy released 
by mass accretion at rate 10~ 8 M Q yr _1 . 



4. DISK MODEL 



We use the minimum- mass model of the protosolar disk (IHayashi et al.lll985l ). with 



1/2 



surface mass density X = 1700 g cm -2 (r/AU) -3 / 2 and temperature T = 280 K (r/AU)" 
at distance r measured in astronomical units from a Solar-mass star. Vertical hydrostatic 
balance yields a volume density p = (Y>/\/2irH) exp(— z 2 /2H 2 ) declining with height z as 
a Gaussian, with a characteristic scale H = c s /Q that is the ratio of the isothermal sound 
speed c s to the orbital frequency Q. 
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5. DEAD ZONE 

Magnetic stresses extract angular momentum from the disk only if the gas is sufficiently 
ionized to couple to the magnetic fields. Upper bounds on the allowed Ohmic resistivity r\ 
come from considering the dissipation scale. The minimum time for regenerating the field is 
the orbital period, while in the case of global magnetic fields, the characteristic size of the 
variations in the field is the distance to the star. The orbital shear acting on a weak radial 
seed field can produce large-scale toroidal fields faster than the fields diffuse away if 



2 

> 10 (5) 



where vk is the Keplerian orbital speed, Q the orbital frequenc y and power-law radial vari- 



ation is assumed for the magnetic fields (ITurner fc Sand 120081 ) . Using a similar approach 



MRI turbulence regenerates tangled magnetic fields if the resistivity is unable to erase mag- 
netic variations over the MRI fastest-growing vertical wavelength 2ttva z /^, where va z is the 
vertical component of the Alfven speed. The criterion for MRI turbulence 



V Az 



depends on the volume-averaged squared vertical Alfven spe ed, whether the background 



magnetic field is vertical, toroidal, or zero (ISano fc Stondl2002l ). Of eqs. [5] and El the turbu 



lence condition is the more restrictive since the MRI wavelength typically is much less than 
the radius. The two criteria divide the disk into three unequal zones. In the active zone 
both conditions hold and strong magnetic activity is expected. In the undead zone the shear 
condition eq. [5] holds and the turbulence condition eq. [6] does not. Magnetic fields can then 
be generated by shear but not by MRI turbulence. In the dead zone, neither condition is 
satisfied and the gas is almost completely decoupled from the magnetic fields. 

We determine the spatial extent of the magnetic activity in the minimum-mass protoso- 
lar disk by finding the locations where eqs. [5] and [6] are marginally satisfied. The ionization- 
recombination reaction network is solved to determine the resistivity on a grid of 60 x 250 
points spaced logarithmically in radius and linearly in z/H, from r = 0.1 AU to 100 AU and 
from z = to 5H. The Alfven speed in eq. [6] is calculated assuming the vertical magnetic 
fields at each radius have a pressure 1000 times less than the midplane gas p ressure. Similar 



press ure ratios are found in ideal-MHD stratified shearing-box calculations (jMiller fc Stone 



2000). 



We first compute the sizes of the undead and dead zones in a fiducial model disk includ- 
ing ionization by stellar X-rays and long-lived radionuclides, and recombination on grains. 
For this calculation only, the chemical network is solved on a finer grid of 1000 x 1000 points 
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Fig. 2. — Dead zone (red) and undead zone (blue) in the minimum-mass protosolar disk with 
well-mixed l-/im dust grains at a mass fraction 1%. The ionization is due to stellar X-rays and 
long-lived radionuclides, and the vertical component of the magnetic field has pressure 0.1% of 
the midplane gas pressure. The vertical mass column exceeds 1, 10 and 100 g cm~ 2 in the light, 
medium and dark gray shaded regions, respectively. Note that the vertical scale is expanded by 
a factor four. 

to more precisely determine the shape of the dead zone boundary. The results appear in 
figure [21 MRI turbulence reaches the midplane only outside 15 AU, while large-scale mag- 
netic fields couple to the orbital shear at the midplane outside 5 AU. The boundary of the 
turbulent region approximately follows the contour of vertical column 10 g cm -2 , close to the 
X-ray absorption depth of 8 g cm -2 . Beyond 12 AU, turbulence occurs at greater columns 
because the lower densities mean slower recombination. 

Results for eight different external ionization scenarios are shown in figure [3J All the 
scenarios include stellar X-ray ionization. The first has no other external sources of ion- 
ization, and corresponds to figure [2j The second includes ordinary interstellar cosmic rays. 
The third and fourth have the cosmic rays enhanced by factors 10 3 and 10 5 . The fifth and 
sixth instead have energetic particles from the disk corona, one with the scaling indicated 
in section 12.41 and the other with an ionization rate 100 times greater, corresponding to a 



- 12 - 



higher mass accretion rate of 10 -6 M yr -1 . The seventh and eighth instead have protosolar 
energetic particles, one with the rate from section 12.51 and the other with a rate 100 times 
greater. All the calculations in figure [3J also include ionization by long-lived radionuclides, 
and recombination on grains. 

The dependence on the radionuclide ionization rates and on the presence of grains is 
shown in the remaining figures. Short-lived radionuclides replace the long-lived radionuclides 
in figure HI The dust is removed for figures through [7J which include no radionuclides 
(figure E]), long-lived radionuclides with an abundance reduced by a factor 10 4 below the 
dusty cases (figure [6]) and short-lived radionuclides with the same abundance reduction 
factor (figure [7]). 

The results with and without dust are quite different. With dust, the active zone fills 
the top 10 to 20 g cm -2 if the external ionization comes from stellar X-rays (figure [3]) • The 
active zone is little changed by ordinary cosmic rays or by energetic particles from the disk 
corona, since these are at most comparable to the X-rays near the bottom of the active zone. 
The active column increases to 20 to 100 g cm -2 if the cosmic ray flux is boosted a thousand 
times, and to 100 g cm -2 if the disk is exposed to the protosolar energetic particles inferred 
from X-ray flares in young stars. The column of material lying in the undead zone is three to 
ten times that in the active zone, while any remaining column is dead. Including short-lived 
radionuclides moves the dead zone boundary slightly deeper inside the disk while having 
little effect on the outer boundary of the undead zone (figure H]). 

With no dust, the active layer fills the outermost 100 g cm -2 if the ionization is due to 
stellar X-rays, increasing to 100-200 g cm -2 if ordinary cosmic rays are added, and reaching 
400-800 g cm" 2 under an enhanced flux of cosmic rays, disk energetic particles, or protosolar 
energetic particles (figure E]) • The undead column is two to ten times that in the active layer. 

Including a small amount of radioactivity leaves the active layer unchanged because the 
radionuclide ionization rate is much less than that due to X-rays at the active layer boundary. 
The extra ionization can be important however in the regions shielded from the energetic 
protons and photons, yielding a bigger undead zone at the expense of the dead zone. The 
differences are slight with 0.01% of the long-lived radionuclides present in the gas (figure E]). 
By contrast, the whole dead zone becomes undead if 0.01% of the short-lived radionuclides is 
included (figure [7]). Reducing the fraction to 0.001% leads to the appearance of a small dead 
zone just outside 0.1 AU. The threshold internal ionization rate for converting the entire 
dead zone into an undead zone is thus about 3 x 10~ 5 (sr ~ 10~ 23 s -1 . 

Cosmic rays from a nearby supernova can increase the active column by up to a factor 
eight over the interstellar cosmic ray case if the base of the active layer is defined solely by 
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Fig. 3. — Column depth to the top of the undead zone (short dashes) and dead zone (long dashes) 
in the minimum-mass protosolar disk with well-mixed 1 /im dust grains at a mass fraction 1%. 
Solid lines show the midplane column. All eight ionization scenarios include stellar X-rays and 
long-lived radionuclides. The scenarios are (left to right and top to bottom): (1) no additional 
ionization, (2) with interstellar cosmic rays, (3) with cosmic ray flux enhanced 1000 times due 
to a nearby supernova, (4) with cosmic ray flux enhanced a further two orders of magnitude, 
(5) with energetic particles from the disk corona, (6) with the disk particles enhanced 100 times, 
(7) with energetic particles from the young Sun, and (8) with the protosolar particles enhanced 
100 times. 



the cosmic ray ionization rate (IFatuzzo et al.l 120061 ) . When the stellar X-rays and the rise in 
recombination rate with depth are included, the higher cosmic ray flux increases the active 
column at 1 AU only from 16 to 41 g cm -2 with dust (figure [3]) and from 210 to 620 without 
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Fig. 4. — Column depth to the top of the undead zone (short dashes) and dead zone (long 
dashes) in the minimum-mass protosolar disk with dust grains and short-lived rather than long- 
lived radionuclides. The eight ionization scenarios are otherwise identical to figure 



(figure [5]), factors of about three. The thousand-fold enhancement in the cosmic ray flux has 
less effect than removing the dust grains. 



6. CONCLUSIONS 

We calculated the locations of the undead and dead zones in the minimum-mass protoso- 
lar disk, under ionization scenarios including stellar X-rays, long- or short-lived radionuclide 
decay, and energetic protons arriving from the general interstellar medium, from a nearby 



- 15 - 




I- i i i i I i i i i I i i i i_J l i i i i I i i i i__L 

-10 1 2-10 1 

log Radius (AU) 



Fig. 5. — Column depth to the top of the undead zone (short dashes) and dead zone (long 
dashes) in the minimum-mass protosolar disk without dust. No radionuclide decay is included. 
The eight ionization scenarios are otherwise identical to figure 



supernova explosion, from the disk corona, or from the corona of the young star. 



The disk contains an undead zone in all the scenarios. All have a dead zone except 
those with the dust removed and a small fraction of radioactive 26 Al remaining in the gas. 
Lacking short-lived radionuclides, even the strongest external ionization process considered, 
the stellar energetic particles, fails to eliminate the dead zone. The thickness of the mag- 
netically coupled layers can increase over the local chemical equilibrium values we have 
computed, however, if ionized material mixes to t he midplane before recombining. Slow 
recombination requires that dust is largely abse nt (llnutsuka fe Sand 120051 ; lllgner fc Nelson 
2006bl : iTurner et al1l2007l : lllgner fc Nelson! hoosl ) . 
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Fig. 6. — Column depth to the top of the undead zone (short dashes) and dead zone (long 
dashes) in the minimum-mass protosolar disk without dust, assuming that 0.01% of the long- 
lived radioactive 40 K remains in the gas phase. The eight ionization scenarios are otherwise 
identical to figure 

Placing these results in the context of existing work, we see that protostellar disks with 
no dead zone can be constructed by 

1. Raising the temperature of the high-column material sufficiently for thermal ionization. 
The temperatures required are likely to be reached through either stellar illumination 
or accretion heating only within a few tenths of an AU of the star, at average mass 
flow rates. At the highest rates inferred for episodic accretion, thermal ionization of 
the interior may extend to a few AU. The central region could alternate between a cold 




I- i i i i I i i i i I i i i i_J l i i i i I i i i i__L 

-10 1 2-10 1 

log Radius (AU) 



Fig. 7. — Column depth to the top of the undead zone (short dashes) in the minimum-mass 
protosolar disk without dust, assuming that 0.01% of the short-lived radioactive 26 AI is found in 
the gas phase. The eight ionization scenarios are otherwise identical to figure [3] The dead zone 
is eliminated under these conditions, but the undead zone remains. 



state with a dead zone and a low accretion rate, and a hot, ther mally-ionized state 
with a higher accretion rate (lArmitage et al.ll200ll ; IZhu et al.ll2009l ). 



2. Redu cing the surface density below the minimum-mass protosolar model (iFromang et al. 



20021 ). Forming planets under these conditions may be problematic. Lower surface den- 



sities could occur however during the dis persal of the disk, for example as X-ray ionized 
material is eroded from the inner edge (jChiang &: Murray- Clayl 120071 ). 



3. Removing the small dust grains while retaining a fraction of the radionuclides in the 
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gas. The removal could occur through settling in the dead zone (ICieslal 120071 ) or 
incorporation into planetesimals. Making grains unimportant compared with gas-phase 
recombination req uires reducing the dust abundance by several orders of magnitude 



(jSano et al.ll2000l ). a difficult feat during planet formation, wh en solid bodies often 



collid e fast enough to shatter and produce secondary particles (e.g. lWetherill fc Stewart 



1993|). 



Each of these three possibilities appears feasible only in limited circumstances. A fourth 
opti on involves ionization by the electrons carrying currents in the weakly-ionized disk inte- 



rior ( jlnutsuka fc Sanoll2005l ). Further study of the microphysics of this process is desirable. 
We conclude that the dead and undead zones are robust consequences of cold, dusty gas 
without strong internal sources of ionization, having mass columns exceeding 1000 g cm -2 . 
Weaker turbulence at the midplane and reduced or fluctuating mass flows are likely under 
these conditions, and should be considered further in planet formation models. 

We thank Eliot Quataert and the staff of the Theoretical Astrophysics Center at the 
University of California, Berkeley for their hospitality and comments, and Eric Feigelson, 
Hugh Hudson and Mark Wiedenbeck for useful discussions. The work was carried out in 
part at the Jet Propulsion Laboratory, California Institute of Technology, with support from 
the NASA Solar Systems Origins Program. Copyright 2009. All rights reserved. 
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